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ON ESTIMATING THE VENUS SPIN VECTOR 
Dr. Peter D. Argentiero 


The Planetary Explorer (PE) mission is designed to land eight probes on 
Venus during the 1977 Venus opportunity. All these probes will be 
equipped with devices for returning range rate data to the Earth. These 
devices are expected to survive impact and to continue transmitting data 
from the planet’s surface. The Venus spin vector as well as the positions of 
the probes on the surface are observable in these data. Thus the correct 
processing of these data should provide at least some improvement in our 
knowledge of these parameters. This prospect is an attractive one because 
of the considerable interest in the spin vector of Venus and because the 
data necessary for the estimation is obtained as a by-product of the PE 
mission and thus its acquisition imposes no further constraints or com- 
promises on mission planning. 

The purpose of this paper is to indicate the improvement in spin vector and 
probe position estimates one may reasonably expect from the processing 
of such data. This was done by duplicating the ensemble calculations as- 
sociated with a weighted least squares with a priori estimation technique 
applied to range rate data that were assumed to be unbiased and uncor- 
related. The weighting matrix was assumed to be the inverse of the co- 
variance matrix of the noise on the data. Attention is focused primarily 
on the spin vector estimation. 


PREVIOUS ESTIMATES OF VENUS SPIN VECTOR 

Several radar determinations of the Venus spin vector have been reported 
in the literature. The technique involves the transmission of a cw signal to 
the planet and the analysis of the power spectrum of the return signal. If 
the transmitted frequency corrected for relative motion between Venus 
and the tracking station is used as the zero point of the spectrum, then 
every other frequency can be interpreted as a doppler shift, and the in- 
tensity associated with that frequency can be related to a line of constant 
radial velocity on the planet. These lines are easy to describe. A plane 
going through the tracking station and parallel to the spin axis of Venus 
and intersecting the planet’s surface describes a circle. That segment of the 
circle visible to the tracking station is a member of the family of lines of 
constant radial velocity. A frequency at which a peak in the spectra 
occurs corresponds to a region on Venus that is rougher than adjacent 
regions. Thus a cw spectrum is a type of map of the radar brightness of 
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the surface of the planet. By obtaining such spectra at different times it 
should be possible to track specific surface features and from such infor- 
mation infer the spin vector of the planet under investigation. Carpenter 
(Ref. 1) reproduces several of these radar spectra and identifies several 
surface features as peaks in the spectra. The mathematics of how one 
might use such information to estimate a spin vector is given in some detail 
by Shapiro (Ref. 2). Basically a weighted least-squares estimation proce- 
dure is used with the surface features associated with peaks in the spectra 
treated as point sources of range rate data. This method has been applied 
by Carpenter (Refs. 3 and 4), Smith (Ref. 5), Goldstein (Ref. 6), and 
Shapiro (Ref. 2). Their various estimates of the rotation period of Venus 
are quite close to the so-called synodic resonance period of 243.16 days. 

If Venus were to have this period, it would rotate backward relative to the 
Earth four times between each inferior conjunction and thus present the 
same side of the Earth at each inferior conjunction. The uncertainties 
associated with present estimates are such that the hypothesis that Venus 
is in synodic resonance is still unsettled. 

The radar determinations of the Venus spin vector are statistically incom- 
patible (Table 1). Clearly the standard deviations give a too optimistic 
representation of the accuracy of the radar determination technique. The 
reason is that ensemble calculations giving standard deviations of the least- 
squares estimation process do not reflect the impact of modeling errors. 

The most obvious modeling error associated with the technique is that of 
treating the source of a peak in a frequency spectrum as a point. Any 
surface feature on Venus capable of causing a significant perturbation in 
the return spectrum of a radar scan must be quite extensive. A more subtle 
though possibly less important modeling error is the assumption that the 
return rate data are uncorrelated. Corrections must be made to the doppler. 
data before they can be used for least-squares estimation. These correc- 
tions are responsible for a portion of the noise on the data. If this portion 
is significant, then it is improper to model the noise on the data as an un- 
correlated or “white” random process. The correct procedure would then 
be to solve for a bias and perhaps a scale factor error in the data. If this 
were done the ensemble calculation would yield a more accurate though 
quantitatively less impressive figure for the standard deviation of the Venus 
spin vector estimate. Another possibility is that the cause of the peaks on 
the return spectra may be a complex combination of phenomena rather 
than a fixed surface feature. 

It is not easy to obtain reasonable estimates of the quality of the radar 
determination technique. Carpenter (Ref. 1) suggests that the formal 
standard deviation numbers associated with radar determinations should 
be increased several times. This is quite vague, but it would be imprudent 
to be more precise. 
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TABLE I— Recent Radar Determinations and Associated Standard 
Deviations of Venus Spin Vector 


Determination 

Rotation period, 
days 

Right ascension, 
deg 

Declination, 

deg 

Carpenter (Ref. 1) 

242.98 ± 0.04 

94 ±3 

-71.5 + 1 

Shapiro (Ref. 7) 

243.09 ±0.18 

84.7 + 1.8 

-65.8 ±1.2 

Dyce, Pettengill, 
and Shapiro (Ref. 8) 

244.3 ±2 

90.9 ± 1 

66.4 ± 1 

Goldstein (Ref. 9) 

242.6 ± 0.6 

98 +5 

-69 ±2 


SPIN VECTOR ESTIMATION WITH PE PROBES 

The PE mission will consist of two separate launchings of multiprobe space- 
craft. Each launch will land a main probe equipped with a two-way doppler 
device and three miniprobes, each equipped with a one-way doppler device. 
The geometric distribution of the probes is given by Figure 1. The succeed- 
ing analysis rests on the following modeling assumptions: 

(1 ) For the duration of their transmission of doppler data, the 
probes remain stationary relative to the Venus surface. 

(2) For this same duration the Venus spin vector is constant. 

(3) The noise on the range rate data transmitted by the probes is 
stationary and uncorrelated. 

(4) Unbiased estimates of the Venus spin vector and the locations 
and effective radii of the probes are available. These estimates 
are statistically uncorrelated and their standard deviations are 
known. 

Assumptions (1) and (2) should disturb no one. Assumption (3) is some- 
what more troublesome. We defer discussion of this assumption to a later 
stage of the analysis. Concerning assumption (4), the estimate of the Venus 
spin vector and its associated uncertainty would presumably be borrowed 
from a radar determination. The probe positions and associated uncertain- 
ties would be obtained from extensive postflight analysis of all relevant 
data gathered during the PE mission. This would include accelerometer, 
temperature and pressure measurements, and range rate data. 







SCIENCE ACCOMPLISHMENTS, 1971 


91 


NORTH POLE 



(a) First launch. 
NORTH POLE 



(b) Second launch. 

Figure 1 -Position of probes on Venus surface. 
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The conventional way to process the doppler data from the PE probes is to 
form tire usual weighted least squares with a priori loss function and choose 
the spin vector and probe positions that minimize this loss function. To be 
specific, let X be a vector of dimension 27 whose elements are estimates of 
the spin vector and positions of each of the eight probes. Let the vector Y 
be the range rate measurements, arranged in some sequence, which would be 
obtained in the absence of noise if X contained the true values of the param- 
eters in question. Y is a known function of X symbolized by 

Y =KX) (1) 

Let Q be the covariance matrix of the noise on the observations and let P be 
the covariance matrix of an a priori estimate X' of the parameters. Then the 
loss function is 

L(X) = [Y -/(X )] T Q 1 [Y -/(X)] + (X' - X) 5 P- 1 (X' - X) (2) 

where Y is the vector of measured values of range rate data obtained from 
the probes. Notice that by assumptions (3) and (4) the matrices Q and P 
are diagonal, a fact that greatly simplifies the following computations. The 
least-squares estimate is defined as that value of X that minimizes the loss 
function of equation (2). Because the data necessary to implement this 
procedure are not yet available, interest is focused on just the statistical 
properties of this estimation procedure. The covariance matrix of the least- 
squares estimator may be obtained if one more assumption is imposed. It 
must be supposed that the function of equation (1) can be accurately repre- 
sented as a first-order Taylor series about X where X is the least -squares esti- 
mate. Thus we assume that equation (1) can be written as 

Y -/(X) = A (X - X) (3) 

The symbol A represents the variational matrix. If N is the total number of 
range rate measurements, then A is an/V- by 27 -dimensional matrix. The 
element in the fth row and /th column of A is the partial derivative of the 
rth component of Y with respect to the /th component of X. It is relatively 
easy to obtain an analytical expression for A. The details are found in 
appendix 1 of Reference 10. If equation (3) is valid, then the covariance 
matrix of the least-squares estimator is given by 

cov (X) = (A t Q 1 A + P -1 )- 1 
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Equation (4) provides a mode for the performance of a parametric study of 
the accuracy attainable in the estimation of the Venus spin vector. Interest 
was focused primarily on the variation of this accuracy with respect to varia- 
tions in the following parameters. 

(1) Length of time the probes survive on the surface 

(2) Quality of a priori information 

(3) Size of the noise on the data. 

This paper reports on the results of such a parametric study. As usual, nomi- 
nal values were established for all relevant parameters. Certain parameters 
were then systematically varied with the other parameters fixed at their 
nominal values. Equation (4) is then used to obtain the corresponding 
standard deviations of the spin vector and probe locations. There is a 
certain amount of arbitrariness involved in the selection of nominal values. 

It is not easy, for instance, to decide what are reasonable values for the 
a priori uncertainties of the spin vector at some time several years in the 
future when the PE mission is to be executed. Also the time at which the 
probes land is a factor because this affects tracking sight geometry. This 
too is somewhat arbitrary although one would of course choose a time 
during which dual coverage from the DSN is possible. The date chosen 
for the landing is May 19, 1977. Dual coverage from Goldstone and 
Madrid is assumed. The nominal value of the spin vector as obtained from 
reference 1 is 242.982 days for the period and 94.1° and -71.4°, respectively, 
for the right ascension and declination of the spin vector (equator of 1950). 
The longitude and latitude of each probe is given in target coordinates in 
reference 1 1 . The data acquisition rate is assumed to be 1 min -1 . It is also 
assumed that all biases are estimated in the least-squares procedure. 

A summary of a priori standard deviations is provided in Table II. 


TABLE II— A Priori Standard Deviations 


Parameter 

Standard deviation 

Right ascension of main probes 

0.115° 

Declination of main probes 

0.115° 

Right ascension of miniprobes 

0.165° 

Declination of miniprobes 

0.165° 

Range rate from main probes 

5 mm s' 1 

Range rate from miniprobes 

5 cm s' 1 

Period of Venus rotation 

1 day 

Right ascension of Venus spin vector 

10° 

Declination of Venus spin vector 

5° 

Effective radius of Venus at each probe 

5 km 
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The usual standard deviation figure for two-way range rate and for a one- 
per -minute sampling rate is 1 mm s' 1 . See for instance Blackshear and 
Williams (Ref. 12) who use this figure for two-way range rate data in an 
error study similar to this one but for the Viking project. The larger 
figure of 5 mm s' 1 for the two-way range rate data on the main probes 
was used to compensate for the fact that the data were assumed to have no 
time correlation. Because the various corrections that must be made to 
the data in fact do tend to introduce time correlations, this assumption is 
a questionable one. Its legitimacy can only be defended if a bias on the 
data is also solved for in the estimation procedure. In this error study, no 
provisions were made for bias estimation. Hence it was felt that to use the 
usual standard deviation figure for the noise on the two-way range rate data 
would lead to an excessively optimistic result. The use of a larger standard 
deviation number was an effort to compensate for this optimism. Because 
the presence of biases on one-way doppler data is a much more serious 
problem, the standard deviation of the noise on the miniprobe data was set 
at 10 times the corresponding figure for the noise on the main probe data. 

Figures 2,3, and 4 display the possibilities of estimating the period, right 
ascension, and declination of the Venus spin vector as functions of the 
length of time the probe range rate devices survive on the surface. Because 
it is not certain that all range rate devices will survive parachute openings 
or impact on the surface, the figures also demonstrate the deterioration in 
the estimation procedure if a miniprobe per launch fails to function and 
also if a miniprobe and a main probe on each launch fail to function. A 
salient feature of the figures is that most of the improvement in the esti- 
mates occurs within the first few minutes after impact. This fact suggests 
that the primary reason for the feasibility of this estimation procedure is 
the coupling provided by equation (4) between the spin vector estimate 
and the relatively small standard deviations of the a priori estimates of 
probe locations. If this were true, one would expect the quality of the 
estimation procedure to be far more sensitive to changes in the standard 
deviations of the a priori estimates of probe locations than to changes 
in the standard deviations of the noise on the data. This appears to be the 
case. For example, with regard to the nominal case with all probes lasting 
30 min, the standard deviation of the rotation period is 0.28 days. When 
the standard deviations of the noise on the data are doubled, the number 
becomes 0.31 days. With the standard deviations of the noise returned to 
nominal values but with the standard deviations of the a priori probe loca- 
tion estimates doubled, the standard deviation of the estimate of the period 
rises to 0.4 days. 
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Figure 2-Standard deviation of period as a function of time. A priori 
standard deviation of period is 1 Earth day. 
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Figure 3-Standard deviation of right ascension as a function of time. 
A priori standard deviation is 10°. 
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Figure 4 -Standard deviation of declination as a function of time. 
A priori standard deviation is 5°. 



98 


GODDARD SPACE FLIGHT CENTER 


It is apparent from a comparison of Figures 3 and 4 that the right ascension 
of the Venus spin vector is far more observable in this experiment than the 
declination. A glance at Table I reveals that the situation is precisely the 
opposite with regard to the radar determination procedure. It is not 
obvious why this should be so. But in this sense, at least, the two estima- 
tion procedures should neatly complement each other. 

The ability of this estimation procedure to improve knowledge of probe 
positions is somewhat less impressive. After 30 min of tracking, knowledge 
of the positions of the main probes is improved by approximately 15 per- 
cent. Knowledge of the positions of the miniprobes after the same period 
of tracking is improved by approximately 10 percent. After an hour of track- 
ing these percentage improvements are, respectively, 18 and 15 percent. 


CONCLUSIONS 

This paper has demonstrated the feasibility of utilizing range rate data 
generated by PE probes to estimate the Venus spin vector. The estimation 
procedure is viable even if each probe transmits doppler data for just a 
few minutes after impact. The standard deviation figures associated with 
this estimation procedure are not dependent on questionable modeling 
assumptions for their validity. Consequently they are a true measure of 
the estimate’s accuracy. 

The right ascension of the spin vector is far more observable in the data 
than the declination. Fortuitously, the radar determination technique 
estimates the declination better than the right ascension. Thus the two 
estimation techniques complement each other. 

The improvements in probe position uncertainty are somewhat marginal. 

Finally, this proposed experiment uses data that is generated essentially 
as a by-product of the PE mission and hence its implementation poses no 
additional constraints on PE mission planning. 
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